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Non-compact QED3 coupled to a four-fermi interaction 
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We present preliminary numerical results for the three dimensional non-compact QED (QED3) with a weak 
four-fermion term in the lattice action. Approaches based on Schwinger-Dyson studies, arguments based on 
thermodynamic inequalities and numerical simulations lead to estimates of the critical number of fermion flavors 
(below which chiral symmetry is broken) ranging from Nf^ = 1 to Nfc = 4. The weak four-fermion coupling 
provides the framework for an improved algorithm, which allows us to simulate the chiral limit of massless fermions 
and expose delicate effects. 



Three dimensional QED is an interesting and 
challenging field theory which serves as a labo- 
ratory to study dynamical mass generation. It 
is confining and it is believed to exhibit chiral 
symmetry breaking when the number of fermion 
flavors Nf is smaller than a critical value Nfc- 
The model is super-renormalizable, because the 
gauge coupling has dimension i and so provides 
an intrinsic scale of the theory. Interest in this 
theory has recently been revived by suggestions 
that QED3 may be an effective theory for the 
underdoped and non-superconducting region of 
the phase diagram of high- Tc superconducting 
cuprate compounds 0. Preliminary numerical 
results of QED3 with Fermi and gap anisotropics 
were presented at this conference 

The principal analytical approach to study dy- 
namical symmetry breaking in QED3 is the self- 
consistent solution of truncated Schwinger-Dyson 
equations for the fermion propagator. Most re- 
cent attempts yielded Nfc ~ 4 0. An argu- 
ment based on a thermodynamic inequality yields 
the prediction Nfc < | j^- Recent attempts to 
decide the issue by numerical simulations placed 
an upper bound on the dimensionless condensate 
/3^{ipip) of 0(10-4) for TV/ = 2 [S]. It was also 
shown that for ^/ = 1 the P'^iip'ip) is of 0(10-3) 
but it appeared difflcult to determine whether 
Nf = 4 lies above or below Nfc [HI- 

Our purpose is to develop firmer control over 



the computer simulation studies of chiral symme- 
try breaking by studying the theory with massless 
fermions. Dealing with almost massless fermions, 
as must be done in the search of continuum and 
chiral limits of lattice models is very difficult with 
the standard algorithms. However, it has been 
shown that when one introduces in {3+l)d QED a 
four-fermion interaction in a semi-bosonized fash- 
ion the performance of the algorithm improves 
dramatically [Jj. The auxiliary field a{x) acts 
as a site dependent dynamical mass term in the 
fermion propagator. In the broken phase, the 
Dirac operator is now non-singular for fermions 
with zero bare mass and its inversion is accurate 
and very fast. 

We have chosen to add to the QED3 lagrangian 
a discrete Z2 chiral invariant four-fermion term, 
which is preferable over four-fermion terms with 
continuous chiral symmetry because the latter are 
not as efficiently simulated due to massless modes 
in the strongly cut-off theory. The lattice action 
using staggered lattice fermion fields Xi X is given 
by the formulas below: 
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Figure 1. Order parameter for Nf — 4 versus (3s 
at fixed /? = 0.2. The lattice size is 12^. 



Figure 2. Dimensionless condensate for Nf = 4 
versus (3s at fixed (3 = 0.2. 
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The indices a;, a;' consist of three integers {xi, X2, xg) 
labelling the lattice sites, where the third direction 
is now considered to be timelike. The symbol {x,x) 
denotes the set of the 8 dual lattice sites x surround- 
ing the direct lattice site x. Since the gauge ac- 
tion is unbounded from above, this action de- 
fines the non-compact formulation of lattice QED. 
The r]fj,(x) are the Kawamoto-Smit staggered fermion 
phases ( — 1)^1"* ^^,^-1^ designed to ensure relativis- 
tic covariance of the Dirac equation in the contin- 
uum limit. Antiperiodic boundary conditions are 
used in the timelike direction and periodic boundary 
conditions in the spatial directions for the fermion 
fields. The phase factors in the fermion bilinear 
are defined by U^ti = wpiiaxij,), where a^i^ is the 
gauge potential. In terms of continuum quantities, 
axfj. = agAf_,{x), /3 = ■4-, /3s = 4- where a is the 
physical lattice spacing, gl is the four-fermion cou- 
pling and is the gauge coupling. The numerical 
results presented here were obtained by simulating 
the QED3 action using a standard Hybrid Molecular 
Dynamics algorithm. 

The vacuum expectation value of the auxiliary field 
a is the chiral order parameter. In the strong gauge 
coupling limit /3 — > 0, it is known rigorously that chi- 
ral symmetry is broken Therefore, as (3 increases, 
for Nf > Nfc there must be a chiral symmetry restor- 



ing phase transition for some finite (3c. A mean field 
form fitting the data would indicate that our data 
are compatible with the hypothesis that the contin- 
uum limit of the theory at the critical point describes 
a free field. Finding evidence for an interacting con- 
tinuum model would be very puzzling indeed. For 
Nf < Nfc, since we believe the order parameter is 
exponentially small in the continuum limit, the relic 
of the transition may persist as a crossover between 
weak and strong coupling behavior with a tail in i/>-(/; 
extending to weak gauge coupling. 

First, we performed simulations with Nf = 4. 
In Fig. 0we plot {|cr|) versus (3s at fixed (3 = 0.2. 
The simulations were performed on 12'' lattices. We 
should note that in the zero gauge coupling limit the 
Nf = 4 three-dimensional four-fermion model un- 
dergoes a phase transition at /3s = 0.86(1) 9 . In 
the weak four-fermi coupling limit all the dynam- 
ics are expected to reside in the gauge and fermion 
fields and the extra four-fermion term just provides 
the framework for an improved algorithm. This be- 
comes clearer in Fig.|21where we see that the dimen- 
sionless condensate /3'^('i/)^) at /3 = 0.2 in the weak 
four-fermion coupling limit saturates to a nonzero 
value ~ 0.0055. This result is consistent with re- 
sults reported in which were obtained using stan- 
dard chiral extrapolations. In Fig. |21 we show data 
generated on 24^ and 32'^ lattices for (a) versus (3 
at fixed (3s = 4.0. The difference between the 24^ 
and the 32'^ data is small, implying that the 32^ 
data are close to the thermodynamic limit. The 
32^ data can be fit to the standard scaling relation 
(a) = const. (/3c - /3)''""'s. Por 0.12 < /3 < 0.16 we get 
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Figure 3. Order parameter for Nj = A versus /3 
at fixed = 4.0. 



Figure 4. Order parameter for Nf 
at fixed (3^ = 2.0. 
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/3mag = 0.46(3) which is consistent with mean fieid 
theory. For 0.15 < P < 0.18 we get /3niag = 0.92(11) 
and Pc = 0.211(6). If this value of /3niag persists for 
couplings arbitrarily close to /3c it will become clear 
that the theory has an interacting fixed point. A de- 
tailed finite size scaling analysis is required to clarify 
this issue. However, we should note that this value 
of /9mag is close to its value (/Smag = 1-1(1) at the UV 
fixed point of the (2 + l)d four-fermion model [O]- Al- 
though /3s = 4.0 is deep in the symmetric phase of the 
"pure" four-fermion model, we still have to study in 
detail the interplay of the QED3 and the four-fermion 
fixed points before we reach a conclusion. The fact 
that QED3 is super-renormalizable may imply that 
the four-fermion fixed point has a much larger domain 
of attraction than the QED3 fixed point. Therefore, 
we may need to perform simulations with smaller Ps 
before we are certain that the simulations lie in the 
vicinity of attraction of the QED3 fixed point. 

We also performed simulations with Nf = 8. With 
fixed Pa = 2.0 we have the following results from 
24^ lattices: For P = 0.11,0.13,0.15 we get /3mag = 
0.49(4) which is consistent with mean field theory. 
However, for P = 0.17, 0.24 we get /3mag = 1.15(8) 
and Pc = 0.272(5). This value of P mag is close to 
its value in the (2 + l)d four-fermion model. We are 
currently performing simulations with Ps = 4.0. Al- 
though our results are very preliminary, the data pro- 
vide some evidence that mean field theory scaling is 
valid in a region where the values of {a) are smaller 
than in the mean field theory region of the Ps = 2.0 
case. In the mean field region {a) is equal to the 
fermion mass. It will be interesting to check whether 



this scenario persists at values of P arbitrarilly close to 
the transition or whether a crossover to a non-trivial 
scaling region occurs at smaller values of (a). 
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